Magnetic Excitations in the S= 1/2 Alternating Chain Compound (VO)2P207 
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Magnetic excitations in an array of (VO)2P207 single crystals have been measured using inelas- 
tic neutron scattering. Until now, (VO)2P2C>7 has been thought of as a two-leg antiferromagnetic 
Heisenberg spin ladder with chains running in the a-direction. The present results show unequiv- 
ocally that (VO)2P2 07 is best described as an alternating spin-chain directed along the crystal- 
lographic 6-direction. In addition to the expected magnon with magnetic zone-center energy gap 
A = 3.1 meV, a second excitation is observed at an energy just below 2A. The higher mode may 
be a triplet two-magnon bound state. 



PACS numbers: 75.10.Jm, 75.10.Hk, 75.30.Ds 

Spin gaps are found in many low dimensional inter- 
acting quantum spin systems. Examples of recent in- 
terest include integer-spin Heisenberg antiferromagnetic 
chains JD , spin-Peierls systems Q and spin ladders HQ . 
The insulating magnetic salt (VO) 2 P 2 7 (VOPO) has 
been widely considered to be an excellent realization of 
the two-leg antiferromagnetic Heisenberg spin ladder §. 
In this Letter we report results from inelastic neutron 
scattering experiments using an aligned array of VOPO 
single crystals. Our results show conclusively that the 
magnetic properties of VOPO are not those of a spin lad- 
der. VOPO is instead found to be an alternating Heisen- 
berg antiferromagnetic chain, with weak ferromagnetic 
interchain couplings. Moreover, the magnetic chains in 
VOPO run perpendicular to the supposed ladder direc- 
tion. The measured alternation parameter describing 
VOPO is in the interesting regime in which a novel two- 
magnon bound state may occur J6|. 

The crystal structure of VOPO is nearly orthorhom- 
bic, with a slight monoclinic distortion so that the space 
group is P2i Q. The room temperature lattice param- 
eters are a=7.73A, 6=16. 59A, c=9.58A, and /3=89.98°. 
The magnetic properties of VOPO arise from 5=1/2 
V 4+ ions situated within distorted VOg octahedra. Face- 
sharing pairs of V06 octahedra are stacked in two-leg lad- 
der structures oriented along the a-axis. The ladders are 
separated by large, covalently bonded PO4 complexes. 
The V-V nearest neighbor distance in the a direction 
is approximately 3.86A, one half of a lattice spacing. 
Along the b ("rung") direction the mean V-V distance 
is 3. 2 A within the ladder. The inter-ladder V-V distance 
across the PO4 groups is 5.lA. The structure is illustrated 
schematically in figure 1. 

The susceptibility of VOPO powder was investigated 
by Johnston et al. Q. Although the susceptibility could 
be accurately reproduced by either a spin ladder (with 
J|l « Jj_) or by an alternating chain [^]]|], the expectation 
that the PO4 group would provide a weak superexchange 
path led to a preference for the spin ladder interpreta- 
tion of VOPO. The proposed spin ladder is illustrated 



schematically in figure 1. Motivated by the detailed the- 
oretical predictions for the excitation spectrum of a spin 
ladder ||, Eccleston et al. || used pulsed inelastic neu- 
tron scattering to probe the dynamic magnetic properties 
of VOPO powders. Their measurement of a spin gap of 
3.7 meV was interpreted as further support for the ladder 
model. 

The surprising discovery of a second spin excitation 
near 6 meV in a recent triple-axis neutron scattering ex- 
periment ]lC| ] on VOPO powder is inconsistent with a 
simple spin ladder model. It has long been known that 
measurements of static magnetic properties cannot eas- 
ily distinguish between spin ladder and alternating chain 
models pd[ , and also that measurements on powders can 
be difficult to interpret. For these reasons we undertook 
studies of spin dynamics in VOPO single crystals. 

Single crystals of (VO)2P207 were grown using a 
method described previously fl(i|j . To achieve sufficient 
sample volume for inelastic neutron scattering measure- 
ments, approximately 200 single crystals of typical size 
1x1x0.25 mm 3 were aligned, mounted on thin aluminum 
plates, and assembled into an array. The resulting sam- 
ple had an effective mosaic spread of 8 — 10° FWHM. 
Although not optimal, this proved to be adequate for 
the present purpose. 

Inelastic neutron scattering measurements were carried 
out using triple axis spectrometers at the HFIR reactor, 
Oak Ridge National Laboratory. The different instru- 
mental configurations utilized are summarized in Table 
I. For temperature control the VOPO array was mounted 
either in a closed cycle refrigerator or in a flow cryostat. 
Both (h, k, 0) and (0, k, I) scattering planes were used in 
different experiments. 

Typical constant-Q scans are illustrated in figure 2. 
The left panels show measurements at (1, —2,0) at tem- 
peratures of 10K (lower panel) and 30K (upper panel). 
The 10K scan clearly shows two peaks, and the reduced 
intensity at the higher temperature is indicative of their 
magnetic origin. The right panel shows a 10K scan at 
Q = (0, —2,0), where the peaks are slightly lower in en- 
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ergy and are much more intense than at (1, —2, 0). 

The solid lines through the 10K data are obtained by 
a least-squares fit of the data to a convolution of the full 
instrumental resolution with delta function excitations 
using dispersion relations discussed below. At low tem- 
peratures the peak widths are resolution limited. The 
significant intensity difference between (1,-2,0) and (0,- 
2,0) is seen to be a consequence of resolution combined 
with the large mosaic spread of the sample. Simple gaus- 
sian fits to obtain the peak positions gave the same results 
as full convolutions. The peak positions were determined 
at various wavevectors. 

Figure 3 shows the measured dispersion along the crys- 
tal axis directions indicated. Several features are imme- 
diately apparent. The excitation energies have at most 
a very weak dependence on the Q c component of the 
wavcvector (middle panel). The dependence of energy 
on Q a (in the "ladder" direction) is much weaker than 
on Qb (in the "rung" direction). This implies that the 
magnetic coupling which supports spin waves is indeed 
dominantly one-dimensional, but that the strongest cou- 
pling is along the b direction! Moreover, the dispersion 
shows that the exchange in the "ladder" direction is fer- 
romagnetic in nature. 

The observed spectrum of magnetic excitations in 
VOPO is clearly inconsistent with a two-leg spin ladder. 
A better model must account for the strong dispersion in 
the b direction, the presence of a spin gap, and be com- 
patible with the crystal structure of VOPO. The most 
likely model is that of weakly coupled antifcrromagnetic 
alternating chains running in the b direction. The gap 
is a consequence of the spin dimerization inherent in the 
VOPO crystal structure. 

A correct theory for VOPO must also explain the pres- 
ence of two magnetic modes. Since the high temperature 
limit of the magnetic susceptibility is exactly that ex- 
pected for S = 1/2, g=2 the possibility that the upper 
mode is a low lying single ion excitation can be ruled out. 

One potential explanation for a second mode is split- 
ting of a triplet magnon level by anisotropy. Using the 
three inequivalent exchanges J\ , J 2 , and J a shown in fig- 
ure 1 the Hamiltonian can be written: 



the anisotropy leads to an S z = longitudinal (L) mode 
and a doubly degenerate S z = ±1 transverse (T) mode: 



H — ^ ^ {Ja a ^Lm &l+l,m + 

l,m a=x, y,z 
jaa qa qa , jaa na na T 



(1) 



where the index l(m) runs over spin sites in the a(b) 
direction. The anisotropic exchange is represented by 
diagonal matrices J? a = J„[5 az + e u (S ax + 5 a v)] for v = 
1, 2, a. Here z is the unique axial direction. 

The excitation energies can be calculated using a 
pseudo-boson technique appropriate for dimerized chains 
|Q| . To first order this agrees with perturbation theory 
from the isolated dimer limit [|13[| . In this approximation 



Uq = {eiJi[eiJi - J 2 cos{TrQ b ) + 2J a cos(irQ a )}} 2 
u5 = {Ji(l + ei)/2[J 1 (l + e 1 )/2- 

e 2 J 2 cos(7rQ b ) + 2e a J a cos(7rQ a )]}5. (2) 

The solid lines in figure (3) are the result of a least- 
squares fit of equation (2) to the data. Two equally 
good descriptions of the data are possible depending on 
whether the upper mode (A) or lower mode (B) is chosen 
as the longitudinal branch; the results are summarized in 
Table II. 

In principle, the Q dependence of the mode intensity 
should allow one to distinguish between the L and T 
modes if the model is applicable. In the present exper- 
iment this is complicated by the large effective sample 
mosaic which dominates the observed changes in mode 
intensity as the direction of Q is varied. Nevertheless, it 
can be stated that both modes appear with similar in- 
tensities for all directions measured. This is unexpected 
in the anisotropic exchange model. 

Fits to the observed dispersion require strong V-V ex- 
change through the PO4 groups. A similarly strong ex- 
change mediated by a PO4 complex has recently been 
found in the chemical precursor compound of VOPO, 
VODPO4 • 1/2D 2 (W), which contains magnetically 
isolated S = 1/2 V-V dimers. Relatively strong su- 
perexchange through the PO4 groups is possible because 
of coherent molecular electron orbitals in the complex. 
Beltran-Porter et al. M| anticipated the importance of 
superexchange through PO4 in VOPO and proposed an 
alternating chain in the 6-direction, which is in agreement 
with the present result. 

Although equation (2) gives an excellent account of 
the data, considerable exchange anisotropy is necessary 
to account for the mode splitting. Recent single crystal 
magnetic susceptibility measurements fl6f| were quanti- 
tatively consistent with the previous powder results || 
and found little, if any, evidence for anisotropy. Fur- 
ther, the exchange coupling in the precursor compound 
VODPO4 • 1/2D 2 was found to be isotropic @. These 
facts suggest that one should seek another explanation 
for the second mode. 

In a recent paper Uhrig and Schulz have consid- 
ered the excitations of an alternating S = 1/2 Heisen- 
berg chain with isotropic exchange (equation (1) with 
£1 = e 2 = 1 and J a = 0) in the continuum limit using field 
theoretical methods. The basic excitations are triplet 
magnons with an energy gap A. The magnitude of A de- 
pends on the alternation parameter 5 = ((3 — l)/(/3 + 1) 
where (3 — J\j J%. The limit (5 = corresponds to the 
usual antiferromagnetic Heisenberg chain for which the 
excitations are free spinons resulting in the well known 
continuum excitation spectrum 17 1. For 8 > a gap 
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opens in the spectrum |6[|lq] and the lowest lying mag- 
netic excitation is a well defined S = 1 triplet mode. For 
small S one expects a second gap to the remnant of the 
original continuum at an energy of 2 A. Both the con- 
tinuum pl| and second gap have been reported in 
CuGeC>3, for which 8 ~ 0.05. The theory further predicts 
that for 5 ~ 0.082 a second sharp triplet mode appears 
at energies just below the continuum |^|. This mode is a 
two-magnon bound state. Assuming isotropic exchange 
the alternation parameter describing VOPO is 6 « 0.1, 
in the regime where this mode should be visible. Given 
the evidence in favor of isotropic exchange in VOPO, it 
is plausible that the observed second mode is in fact the 
two-magnon bound state. 

Figure 4 shows an expanded view of the dispersion in 
VOPO along the chain direction near Q = (0,-2,0), 
which corresponds to Q = it in reduced units. The 
solid and dashed lines are the result of fitting the one- 
dimensional lower (I) and upper (u) mode dispersions to 
the expression 



= A /(^J)W(Q) + A 



2 

/ , u 



(3) 



resulting in J = 9.3(1) meV, A; = 3.12(3) meV, and 
A u = 5.75(2) meV. The expression for the lower mode 
in equation (3) agrees with the predictions of the contin- 
uum model H for Q w 7r. The value of A; /J allows a 
determination of S using the results of numerical compu- 
tations j|. Note that J = (J x + J 2 )/2. The value of A u 
is slightly less than 2A; as expected for the two-magnon 
bound state. In contrast, if the second mode is a conse- 
quence of anisotropy this proximity is purely accidental. 
We believe that a two-magnon bound state is the best 
explanation for the upper magnetic mode. 

Several open questions about this interpretation re- 
main. The substantial intensity of the upper mode is 
very unusual for two-magnon scattering. As indicated 
schematically in figure 4 there may also be a higher lying 
continuum of excitations but no clear sign of this has ap- 
peared in our measurements to date. Detailed theoretical 
calculations of the expected neutron scattering response 
would be of great interest. Light scattering measure- 
ments may be sensitive to two-magnon states and should 
also be instructive. 

In the presence of a sufficiently strong external mag- 
netic field the predicted mode splittings of the two pro- 
posed models are very different. Unfortunately such ex- 
periments are extremely difficult with the VOPO samples 
available to date. A suitably large single crystal is nec- 
essary for this as well as for careful investigations of the 
mode intensities as a function of Q and T. Continuing 
efforts are underway to obtain the desired sample. 

We would like to acknowledge useful interactions with 
P. Dai, J. Fernandez-Baca, H. Mook, and J. Zarestky. 
We thank J. Thompson for sharing results of susceptibil- 
ity measurements with us prior to publication. Expert 



technical assistance was provided by Scott Moore and 
Brent Taylor. This work was supported in part by the 
United States Department of Energy under contract DE- 
FG05-96ER45280 at the University of Florida, and by 
Oak Ridge National Laboratory, managed for the U.S. 
D.O.E. by Lockheed Martin Energy Research Corpora- 
tion under contract DE-AC05-96OR22464. 

TABLE I. Spectrometer configurations^ 



Instrument 




Collimation 




Energy 


HB1A 


40' 


- 40' 


-40' - 


68' 


£/=14. 7meV 


HB1 


48' 


- 40' 


-60' - 


120' 


£ F = 13.5meV 


HB3 


48' 


- 40' 


- 60' - 


120' 


E I ,E F =30.5meV 



t Using PG(002) monochromator, analyzer, PG filter. 



TABLE II. Anisotropic exchange parameters 



Parameter 



fit A 



fit B 



J a (meV) 

ta 

Ji (meV) 
£1 

J 2 (meV) 

£2 



-0.27(2) 
0.91(7) 

10.30(14) 
1.16(3) 
8.7(3) 
1.12(2) 



-0.24(2) 
1.13(9) 

12.32(55) 
0.93(5) 

10.18(15) 
0.85(4) 
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Figure Captions 



Figure 1 (color). Schematic depiction of the struc- 
ture and magnetic interactions in VOPO. The spin ladder 
model previously thought to describe VOPO has nearest 
neighbor exchange constants Jy along the a ("ladder") 
direction and J± along the b ("rung") direction. In the 
alternating chain model, nearest neighbor V 4+ ions are 
alternately coupled by constants J\ and J 2 along the b 
(chain) direction. Neighboring spins in adjacent chains 
are coupled by J a - Magnetic coupling in the c direction 
is negligible. 

Figure 2. Typical constant-Q scans at Q = (1,-2,0) 
(left panels) and Q = (0,-2,0) (right panel). The 
solid lines running through the 10K data are least-square 
fits of a convolution of the instrumental resolution with 
delta function excitations obeying the dispersion rela- 
tions of equation (2). The solid line through the 30K 
data is a guide to the eye. Intensities are normalized to a 
fixed number of incident neutrons expressed by monitor 
count units (MCU). Some representative la error bars 
are shown. 



Figure 3. Measured dispersion of magnetic excitations 
in VOPO at T = 10K. When not visible error bars are 
smaller than the size of the plotted symbols. Filled cir- 
cles (open diamonds) are points from the lower (upper) 
energy mode. The solid lines are dispersion curves calcu- 
lated using parameters obtained by fitting equation (2) 
to the observed peak positions. Wavevectors are plotted 
in units corresponding to the VOPO reciprocal lattice. 

Figure 4. Closeup of the measured dispersion in VOPO 
in the chain direction near the antiferromagnctic zone 
center (0,-2,0), (Q = n in reduced units). The lines 
through the upper and lower modes are from a fit of 
equation (3) to the data for wavevectors along the chain 
direction alone. The hypothesized continuum is illus- 
trated schematically. 
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